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I. INTRODUCTION 



Top quark pair production in pp collisions is a sensi- 
tive probe of quantum chromodynamics at high energy. 
At lowest order in the standard model (SM), quark pair 
production is symmetric under charge conjugation. At 
next-to-lcading order (NLO) the interference of processes 
that differ under charge conjugation leads to a small 
forward-backward asymmetry of order 0.06 ± 0.01 in the 
tt rest frame [IHl] . An analogous effect is predicted at 
order in QED and is confirmed in measurements of 
e+e~ M"*"/^" 15 • Study of the NLO QCD asymmetry 
in inclusive jet events is hampered by the difficulty of 
measuring the jet charge. In pair produced top quarks 
with one semi-leptonic decay, the top can be tagged ac- 
cording to the well-measured lepton charge, enabling a 
probe of the NLO QCD effect and a test of charge con- 
jugation symmetry in strong interactions at high energy. 

The CDF and DO experiments have made initial mea- 
surements of the ti asymmetry in pp collisions at y/s — 
1.96 TeV at the Fermilab Tevatron In the CP- 

invariant pp system, the NLO QCD effect appears as a 
charge dependent forward-backward asymmetry of the 
top quark direction with respect to the proton direc- 
tion. Using data samples corresponding to 1.9 fb~^ 
and 0.9 fb~^ respectively, CDF and DO report pos- 
itive asymmetries that are consistent with the QCD 
prediction within large experimental uncertainties. A 
number of theoretical papers suggest interesting new 
physics mechanisms including axigluons, diquarks, new 
weak bosons, and extra-dimensions that can all pro- 
duce forward-backward tt asymmetries [71 [5] . The model 
building must accommodate the observed consistency of 
the tt cross-section and total invariant mass distribution 
with the SM QCD prediction PEU]. 

We report here on a new study of the forward- 
backward asymmetry in pp collisions at ^/s = 1.96 TeV, 
using data corresponding to an integrated luminosity of 
5.3 fb~^ recorded with the CDF II Detector. We study 
events with the lepton-fjets topology, where either the 
t or i has decayed semileptonically. The asymmetries 
are measured in two variables: yh, the rapidity of the 
hadronically decaying top quark, corresponding to the 
top rapidity in the laboratory (lab) frame, and Ay, the 
difference of the rapidities of the leptonic and hadronic 
top systems, which measures the top quark rapidity in 
the tt rest frame. We show that t and t asymmetries are 
consistent with CP conservation, and combine them to 
measure the total asymmetry in the sample. We measure 
the inclusive asymmetries, and the functional dependence 
of the tt frame asymmetry on Ay and on the total invari- 
ant mass of the tt system, M^f. We apply corrections 
for backgrounds, acceptance, and resolution to calculate 
parton level measures of the inclusive asymmetry in both 
the lab and tt rest frames, and in two regions of Ay and 
Mff in the tt frame. 



II. DETECTION, EVENT SELECTION AND 
RECONSTRUCTION 

CDF II is a general purpose, azimuthally and 
forward-backward symmetric magnetic spectrometer 
with calorimetry and muon detectors jll) . Charged parti- 
cle trajectories are measured with a silicon-microstrip de- 
tector backed by a large open-cell drift chamber in a 1.4 T 
solenoidal magnetic field. Electromagnetic and hadronic 
calorimeters located outside the solenoid provide jet and 
missing energy reconstruction. Outside the calorimeter 
are multilayer proportional chambers and plastic scintil- 
lator hodoscopes that provide muon identification in the 
pseudorapidity region | r/ |< 1.0. We use a cylindrical 
coordinate system with origin at the detector center and 
z-axis along the proton direction [T^] . 

This measurement uses tt candidate events in the "lep- 
ton-|-jets" topology, where one top quark decays semilep- 
tonically {t — )■ Ivb) and the other hadronically {t — > qq'b) 
[13 ■ We detect the lepton and four jets from top quark 
decays and quark hadronization, and an inferred neu- 
trino based on the presence of missing energy. The de- 
tector is triggered by a high transverse momentum elec- 
tron(muon) in the central portion of the detector with 
Exipr) > 20 GeV(GeV/c) and |?7| < 1.0. We require four 
or more hadronic jets with Et > 20 GeV and |r/| < 2.0, 
and a large amount of missing transverse energy, > 20 
GeV, consistent with the presence of an undetected neu- 
trino. The jets are reconstructed using a cone algorithm 
with SR = \J bf^P' -I- St]^ < 0.4, and calorimeter signals 
are corrected for detector inefficiencies and for the en- 
ergy scale factor. The SECVTX algorithm [T3] is used 
to find displaced &-decay vertices using the tracks within 
the jet cones, and at least one jet must contain such a 
"6-tag". Jets with 6-tags are restricted to |?7| < 1.0. 

A total of 1260 events pass this selection. The size of 
the non-tt background processes in the lepton+jets+6- 
tag selection is derived in precision measurements of the 
tt production cross-section [91. The estimated back- 
ground in our sample 283.3 ± 91.2 events. The pre- 
dominant backgrounds are from QCD-induced W-|-multi- 
parton events containing either 6-tagged heavy-flavor jets 
or errantly tagged light-flavor jets. These are mod- 
eled using a simulation sample derived from the ALP- 
GEN generator [15] and a data driven technique that 
derives tagging efficiencies, mis-tagging rates and sam- 
ple normalizations from direct measurement. A back- 
ground component from QCD multi-jet events with fake 
leptons and mis-measured is modeled using multi-jet 
events with lepton candidates that are rejected by our 
cuts. Other small backgrounds from electroweak pro- 
cesses {WW, WZ, single-top) are reliably estimated us- 
ing Monte Carlo generators. Further details on the sam- 
ple selection and background modeling can be found in 
Ref. [Sj. 

The reconstruction (reco) of the tt kinematics uses the 
measured momenta of the lepton, the M^t, and the four 
leading jets in the event. The jet-parton assignment and 
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calculation of the tt four- vectors uses a simple x^-based 
fit of the lepton and jet kinematics to the ti hypoth- 
esis, allowing the jet energies to float within their ex- 
pected uncertainties, and applying the constraints that 
Mw = 80.4 GeV/c2, Mt = 172.5 GeV/c^, and 5-tagged 
jets are associated with b-partons. This algorithm is well 
understood in the context of precision top mass measure- 
ments, where the fit is performed without the top mass 
constraint [H] , and other top physics studies that use the 
top mass constraint [TOl . We study the reconstructed top 
quark rapidity and the difference in the reconstructed top 
and anti-top rapidities, from which we derive the forward- 
backward asymmetries in the pp (laboratory) rest frame 
and in the ti rest frame. 

The validity of the analysis is checked at all steps 
by comparison to a standard prediction made using the 
PYTHIA [T7| tt model, the CDF lepton-|-jets-|-6-tag back- 
ground model, and a full simulation of the CDF-II de- 
tector. We use PYTHIA 6.2.16 with CTEQ5L parton dis- 
tribution functions [18] and Mt = 172.5 GeV/c^. The 
background model developed in concert with the preci- 
sion cross-section studies provides good measures of both 
the normalizations and shapes of the non-tt processes [3] ■ 
The veracity of the combined pythia plus background 
model, and in particular, its reliability for the estimation 
of systematic uncertainties, is well verified in many other 
top-physics studies at CDF [3 El UHl [HI EI] . 

Note that because pythia does not include the NLO 
QCD charge asymmetry, the standard pythia predic- 
tion is not the SM prediction for the forward-backward 
asymmetry. Studies with the mc@nlo generator [TO] (see 
Sec. IV B ) predict that the magnitude of reconstructed 



QCD asymmetry in our sample is smaller than the cur- 
rent experimental resolution. Symmetric pythia is thus 
a good approximation for SM tt and provides an unbi- 
ased control sample for many of our studies. We will com- 
pare our measurements to the SM predictions of mc@nlo 
when appropriate. 



III. RAPIDITY VARIABLES AND 
ASYMMETRY DEFINITIONS 



t quark, {t quark, yf^)- If CP is a good symmetry, 
the distributions of y^^ and yf^ are reflections of each 
other, and we can combine both samples, weighting with 
the lepton charge, to use —qyh as the rapidity of the t 
quark in the lab frame, y^^. 

A frame independent measurement is available in the 
rapidity difference of the leptonic and hadronic systems 
^yih = yi ~yh- After multiplication by the lepton charge 
q, this variable measures the difference between the top 
and antitop rapidities: q/^yih = q{yi - yh) = yt - yi = 
Ay. The rapidity difference Ay is independent of the tt 
system longitudinal motion and is simply related to the 
top quark rapidity in the ti rest frame: y'* = ^Ay. Since 
the rapidity preserves the sign of the production angle, 
an asymmetry in Ay is identical to the asymmetry in the 
top quark production angle in the ti rest frame. 

With N as the number of events with a given rapidity, 
we define the total ti frame asymmetry: 



iV(Ay > 0) - A^(Ay < 0) 
Ar(Ay > Q) + N{/^y < 0) 
Njyf > 0) - Njyf < 0) 
iV(yf > 0) + Niyf < 0) 



(1) 



and the total laboratory frame asymmetry, assuming CP 
invariance: 



(2) 



N{-qyk > 0) - N{-qyh < 0) 
N{-qyh > 0) + N{-qyh < 0) 
jV(yPP > 0) - Njyf < 0) 
iV(yPP > 0) + N{yf < 0) 



Since y^ and Ay;/j are identified with either a t or an i 
by the sign of the lepton in the event, they are the pri- 
mary variables for defining the charge dependence of the 
asymmetries and testing for CP invariance. We define 
the forward-backward charge asymmetry in the ti rest 
frame to be: 



In the lepton-|-jets decay topology of the ti pair, there 
is a leptonic decay, t — >■ Wb Ivh, and a hadronic de- 
cay t — >■ Wh — qq'h. The complications of the central 
lepton acceptance and the reconstruction of the neutrino 
from the J^t create a difference in the reconstruction res- 
olution for the two different kinds of decay. In order to 
control effects of this kind, our treatment of top rapidity 
variables maintains the distinction between the leptonic 
and hadronic decay systems, with the t and i assignments 
following in accordance with the lepton charge. 

The most direct measurement of the top direction with 
respect to the beamline is the rapidity of the hadronic top 
system in the lab frame, yh, which has acceptance out to 
|ry| = 2.0 and good directional precision. In events with 
a negative (positive) lepton, yh is the lab rapidity of the 



± ^ A^^(Ay;;, >0)-iV^(Ay;;, <0) 

N±{l^yih>Q)+N±{^yih<Q) 



and in the laboratory frame to be: 



,± ^ iV^(y,. >0)-jV^(y,,<0) 
" 7V±(y„>0)+7V±(y, <0) 

where the superscript refers to the sign of the lepton 
charge q. 
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TABLE I: Summary of rapidity variables and asymmetries, 
definition 

yh rapidity of hadronic top system in lab 

yi rapidity of leptonic top system in lab 

Ayih rapidity difference yi - yu 

Ay tt rapidity difference: yt - yt = q{yi ~ yh) 

top quark rapidity in laboratory frame: —qyn 
top quark rapidity in ti rest frame: |Aj/ 

Ah asymmetry in yh 

Aih asymmetry in Ayih 

A'^^ laboratory frame asymmetry in yf^ (both charges) 
A*'^ ti frame asymmetry in y" (both charges) 



The laboratory and tt frame present trade-offs for the 
asymmetry measurement. The laboratory frame is exper- 
imentally simple: the direction of the three-jet hadronic 
top decay in the detector is well-resolved, with uncer- 
tainty dominated by a Gaussian width Syh ~ 0.034, and 
free from the complications of the neutrino reconstruc- 
tion [20] . The yh distribution is thus the simplest way to 
test for the presence of an asymmetry. However, as the 
laboratory frame includes an uncontrolled longitudinal 
boost from the rest frame of the primary qq interaction, 
the information on the fundamental production asymme- 
try is diluted. 

Because the momentum scale of initial state radia- 
tion is small compared to Mjf, the qq frame is well ap- 
proximated by the tt rest frame. We measure the tt 
frame rapidity in an experimentally robust way using 
the difference of two rapidities in the detector frame, 
= Qiui — Vh)- But the inclusion of yi and the poorly 
resolved neutrino reconstruction degrades the precision: 
the Gaussian part of the tt frame resolution has width 
5Ay ^ 0.100 and significant non-Gaussian tails. The tt 
frame has an advantage in interpretation, but a disad- 
vantage in resolution. 



The frame dependent resolution has to be considered 
against a possible frame dependence in the size of the 
asymmetry. In the case of the QCD charge asymmetry, 
our NLO models (see Table III) suggest that the recon- 
structed asymmetry is reduced by a factor of 0.6 — 0.7 
in the transition from the tt to laboratory frame. This 
roughly balances the resolution difference to give compa- 
rable sensitivities to the inclusive QCD asymmetries in 
the two frames. Asymmetries generated by other pro- 
cesses may produce a different ratio between the two 
frames, possibly with a Ay or M^j dependence, and a 
more precise measurement of the ratio could help to il- 
luminate the underlying physics. We will return to this 
issue in Sec. IVIIIEI 



A summary of the rapidity variable and asymmetry 
definitions used in this paper is given in Table [l] 



IV. PHYSICS MODELS AND EXPECTATIONS 

We briefly describe the theoretical basis for the QCD 
asymmetry at NLO, the calculation of the theoretical 
asymmetry using the MCFM program and use of 
the MC@NLO event generator in creating a simulated 
NLO sample for input to our analysis. We also de- 
scribe a simple chiral color-octet model, executed in mad- 
GRAPH [23^, that we use to understand the response of 
our analysis to a large tt asymmetry. 



A. NLO QCD Theory and MCFM 

The NLO QCD asymmetry arises in the interference 
of qq processes that behave differently under charge con- 
jugation. The gg initial state does not contribute to the 
asymmetry, but does dilute the average value. 

Early, pre-top, treatments of the QCD charge asym- 
metry discussed measurement of generic heavy quarks 
in hadron collisions (24]. More recent treatments have 
focused on the particular case of the top quark at the 
Tevatron and at the LHC [IHS]. 

The asymmetry gets a positive contribution from inter- 
ference of the tree-level and box diagrams, as in the upper 
diagrams in Fig. [l] and a negative contribution from the 
interference of initial and final state radiation in tt -\- jet 
[tij) final states, as in the lower diagrams. The total 
inclusive asymmetry is the sum of these opposing contri- 
butions. An intuitive picture of the first effect is that the 
QCD Coulomb field of an incoming light quark repels the 
t quark to larger rapidities, while attracting the t quark 
to smaller rapidities, thus creating a positive asymmetry 
at large r], as defined by the quark direction |2T]. The 
second effect can be pictured in terms of color-flow: when 
a forward gluon is radiated by the incoming quark, the 
large acceleration of the color charge biases the top quark 
towards the backward direction. 

Predictions for the NLO QCD asymmetry are derived 
using version 5.7 of MCFM with CTEQ6.1(NL0) ^ and 
Mt — 172.5 GeV/c^. The forward-backward asymmetry 
in the tt rest frame is found to be ^" = 0.058 ± 0.009. 
In the laboratory frame the top quark rapidities are 
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FIG. 1: Interfering qq — >■ tt (above) and qq — > ttj (below) amplitudes. 



broadened by the varying boost of the tt system along 
the beamline, and the asymmetry is diluted to A^p — 
0.038 ± 0.006. Our mcfm predictions are in accord with 
other recent calculations [TH3]. These predictions are for 
top quarks as they emerge from the qq collision, before 
any modifications by detector acceptance and resolution. 
We will call this the parton-level. Based on our own stud- 
ies of scale dependence in mcfm and also the studies in 
the references above, we assign a 15% relative uncertainty 
to all NLO MCFM predictions. 

An NLO calculation for inclusive tt production is an 
LO calculation for the production of a + jet final state, 
and thus an LO calculation for the asymmetry in final 
states containing an extra jet. A new NLO calculation 
for ttj production (and thus for the asymmetry) suggests 
that the negative asymmetry in this final state is greatly 
reduced from leading-order [55]. This new result for the 
ttj asymmetry can be incorporated into an analysis of 
the asymmetry for inclusive tt production only within the 
context of a full NNLO calculation of tt production. Such 
calculations are underway but are not complete. Thresh- 
old resummation calculations indicate that the inclusive 
asymmetry at NNLO should not differ greatly from that 
predicted at NLO [TJ l21j. In this paper, we compare 
to the NLO predictions for ti production. We include a 
15% scale dependence uncertainty, but note that there is 
an overall unknown systematic uncertainty on the theo- 
retical prediction pending the completion of the NNLO 
calculation. 



B. NLO QCD Simulation with MCONLO 

We use the event generator mc@nlo to create a sim- 
ulated sample that includes the QCD asymmetry as pre- 
dicted by the standard model at NLO. In addition to 
including the asymmetric processes this generator prop- 
erly estimates the amount of gg, and thus the dilution of 
the asymmetry from these symmetric processes. 

Some naming conventions for the data-to-simulation 
comparison are given m Table [n] AU Monte Carlo (MC) 
generators will have the same conventions: the truth in- 
formation is the parton level; the pure top signal after 
simulation, selection, and reconstruction is the tt level, 
and the full prediction including backgrounds is tt + bkg 
level. The reconstructed lepton-|-jets sample is the data. 
Subtracting the backgrounds from the data yields the 
reconstructed tt signal-level. Correcting the data for ac- 
ceptance and resolution produces a measurement at the 
parton-level. 

TABLE II: Naming conventions for data and simulation sam- 
ples. 



sample 


level 


definition 


comparable to 


data 


data 


reco 1+jets 




data 


signal 


data minus bkg 


tt in data 


data 


parton corrected signal 


ti at creation 


MC 


tt+bkg 


reco tt + bkg 


data 


MC 


tt 


reco ti no bkg 


data signal 


MC 


parton 


truth level 


data parton 



In the near-threshold form of the cross section T] the 
tt frame asymmetry can be seen to increase with the top 
quark production angle and velocity (/?), and these are 
thus key variables for understanding the source of the 
asymmetry. In this analysis, the proxies for these vari- 
ables are the top quark rapidities and the mass Mjf of 
the ti system. Measurements of the rapid ity and mass 
dependence of A*' are described in Sections VI and VII 



The MC@NLO predictions for the asymmetries at var- 
ious levels of simulation are shown in Table IIIII The 
uncertainties include the Monte Carlo statistics and the 
NLO theoretical uncertainty. The parton-level mc@nlo 
asymmetries are consistent with mcfm, as expected. Af- 
ter CDF detector simulation, event selection, and recon- 
struction, the asymmetries in the mc@NLO ti signal are 
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TABLE III: NLO QCD asymmetries in two frames. Uncertainties includes MC statistics and scale dependence. 

model level A^^ yl" ^pp/tI" 

MCFM parton 0.038 ± 0.006 0.058 ± 0.009 0.66 ±0.10 

mcOnlo parton 0.032 ± 0.005 0.052 ± 0.008 0.62 ± 0.09 

MC@NLO ti 0.018 ± 0.005 0.024 ± 0.005 0.75 ± 0.11 

MC@NLO tt+bkg 0.001 ±0.003 0.017 ±0.004 0.06 ± 0.01 




FIG. 2; Rapidity distributions in data compared to predictions. The legend "tt ± bkg" implies totals in those bins are the sum 
of the tt and background components. The asymmetries in the data and the predicted tt signal, background, and combination 
are shown in legends on top right of plots, using the conventions of Table [III 



significantly reduced. In the laboratory frame, the ex- 
pected asymmetry at the reconstructed tt+bkg level is 
consistent with zero. 

We will see in Sec. |v]that the statistical err or o n A^^ 
and in the current dataset is 0.028. Table III shows 



that, even after background subtraction, the central val- 
ues of the expected asymmetries are smaller than the 
experimental resolution. This motivates the continued 
use of symmetric pythia as our default tt model (as dis- 
cussed in Sec.|ll|, but we will also consider the mc@nlo 
predictions in several specific studies. 



C. Generic Color-Octet with madevent 

It is important that we test our measurement proce- 
dures in the regime of the observed asymmetries. We 
have used madgraph and the model of Ref. [5] to create 
asymmetric test samples that can be used as input to our 
analysis A massive axial color-octet G mixes with 

the gluon to give a production cross section including pole 
and interference terms linear in cos(6'*), where 9* is the 
t production angle in the ti rest frame. In these models 
the asymmetry is an explicit function of the production 
angle and (f' , again illustrating the importance of the Ay 



and Mil dependence for understanding the source of the 
asymmetry. 

We tuned the octet mass Mq to put the pole out of 
range and the couplings to give inclusive parton level 
asymmetries in rough agreement with the data, A^^ = 
0.110 and A" — 0.157, while minimizing the effect on 
the ti cross section and Mjf distribution (see Appendix) . 
After MADGRAPH generation, partons are showered with 
PYTHIA and the sample is passed through the complete 
CDF-II detector simulation. We call this sample OctetA. 

A second sample, OctetB, has the same couplings and 
lower Mq, to give inclusive parton level asymmetries 
= 0.205 aiid A" = 0.282, with shght (~ 5%) in- 
creases in the ti cross section and in the high Mti tail. 
We consider that OctetA is a better model for under- 
standing the experimental response, but we will appeal 
to both models in order to span an asymmetry range 
extending beyond the experimental values. 

We emphasize that our use of the Octet models is to 
study sensitivities and systematic effects in the presence 
of large asymmetries, and should not be construed as 
tests of physics hypothesis. More detail on these samples 
can be found in the Appendix. 
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FIG. 3: Distributions of Ayih (top) and yh (bottom) for events with negative leptons (left) and positive leptons (right). 



V. MEASUREMENT OF THE INCLUSIVE 
ASYMMETRIES 

We now turn to the rapidity distributions in the data. 
The inclusive distributions of the Ayih and yh variables 
are shown in Fig. [2] compared to the standard pythia 
tt + background prediction. These distributions contain 
the full sample of both lepton signs and should be sym- 
metric. The legend on the top right shows the asym- 
metries in all components. The data agrees well with 
tt+hlig prediction in both variables, and, in particular, 
the asymmetries are consistent with zero. 

A forward-backward asymmetry becomes apparent 
when the sample is separated by charge. The top row 
of Fig. [3] shows the Ay distributions for events with neg- 
ative leptons (left) and positive lcptons(right). Wc find 
= 0.067 ± 0.040 and A" = -0.048 ± 0.039, where 
the uncertainties are statistical only. With limited sig- 
nificance, the asymmetries are equal in magnitude and 
opposite in sign. 

The bottom plots of Fig. [3] shows the yh distribu- 



tions for events with negative leptons (left) and pos- 
itive leptons (right). An indication of asymmetry is 
also observed in this figure: t quarks are dominant in 
the forward (proton) direction and the t quarks in the 
backward (p) direction. The measured asymmetries are 
A+ = -0.070 ± 0.040 and Aj; = 0.076 ± 0.039, again 
equal and opposite within uncertainties. 



The sign reversal of the asymmetry under interchange 
of the lepton charge (or, in our formalism, under inter- 
change of t and t) is consistent with CP conservation. 
With larger samples and improved precision, the com- 
parison of the charge separated distributions will pro- 
vide a strict test of CP conservation in tt production. If 
we assume CP conservation we can calculate the total 
asymmetry in each frame using Eqs. (1) and (2). The 
distributions of these variables are shown in Fig. |4j The 
asymmetry in the tt frame is A*' = 0.057 ± 0.028, and in 
the laboratory frame is ^pp = 0.073 ± 0.028, where both 
uncertainties are statistical. 
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A. The Parton-Level Asymmetry 

In order to compare our results to theoretical predic- 
tions we must correct the data for backgrounds, for in- 
complete detector acceptance, and for the finite rapidity 
resolution of the reconstruction. 

We derive the signal level tt distributions by sub- 
tracting the expected background from the reconstructed 
data. This correction is most important in the laboratory 
frame, where, as shown on the right in Fig. [4] the back- 
grounds show a significant negative asymmetry originat- 
ing in the W production asymmetry in VF-|-jets events. 

The reliability of the background model is verified in 
the subset of the lepton-|-jets selection that has no b- 
tagged jets. This "anti-tag" sample is background en- 
riched, with S:B ^ 0.3, and is also fully modeled in our 
analysis. The predicted ti and lab frame asymmetries 
in the anti-tag data sample are in excellent agreement 



with observation, as shown in Table IV The absence of 



asymmetry in this background enriched sample, and the 
consistency between prediction and observation, suggest 
that the asymmetry in the 6-tagged sample is correlated 
with the tt signal and not the backgrounds. 



TABLE IV: Asymmetries in the anti-tag sample of the data 
and ti + bkg level prediction. 



selection 



A' 



anti-tag data 0.033 ± 0.018 -0.016 ± 0.018 

anti-tag prediction 0.010 ± 0.007 -0.023 ± 0.007 



Acceptance and resolution corrections are made with 
a simple linear unfolding of the Ay and yf^ distribu- 
tions using the technique described in Ref. [5]. Let the 
binned parton-level rapidity distributions be represented 
by the vector n. The n distribution is modified by the 



acceptance and then by the smearing in the reconstruc- 
tion. These transformations can be expressed as matri- 
ces transforming the distribution vector from the parton 
level to our reconstructed signal: risignai = SArJparton- 

The matrices A and S are derived from pythia sam- 
ples by comparing distributions at the Monte Carlo truth 
level to the same distributions after reconstruction. The 
acceptance matrix A is diagonal. The smearing matrix 
S measures the bin-to-bin migration arising from the fi- 
nite resolution of reconstructing the events in the tt hy- 
pothesis. To measure the parton-level value, we subtract 
backgrounds to recover the signal from the data, and then 
invert the transformation: 



''parton 



(,'T'data 



"^bkg) 



(5) 



With the assumption of the A and S response as com- 
puted with PYTHIA, this technique gives a model inde- 
pendent measure of the parton-level asymmetry. The 
result was found to be robust and the uncertainty mini- 
mized when the distributions are separated into four bins 
with bin edges at 0.0, and \Ay\ = 1.0 or \yf^\ = 0.5 

[3 [Msg. 

The measurement is affected by uncertainties in our 
models for the amount and shape of the backgrounds, 
the amount of initial state and final state radiation (ISR 
and FSR) in pythia, the jet energy scale (JES) of the 
calorimeter, the parton distribution functions (PDF), 
and the color reconnection in the final state. These ad- 
ditional systematic uncertainties are studied by repeat- 
ing the analysis with reasonable variations in the model 
parameters. We also test the result of substituting the 
other LO generators herwig and ALPGEN for pythia 
in the model for the matrix unfold. The effect of these 
model variations on the parton-level asymmetry is small, 
as seen in Table IVl 
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Ay = y^-y. -qy, = yr 

FIG. 5: Four-bin representation of rapidity distributions for all correction levels. Solid histogram is the pythia tt model. 



TABLE V: Systematic uncertainties on parton-level asymme- 
tries in both frames. 



effect 



background magnitude 0.015 0.011 

background shape 0.014 0.007 

ISR/FSR 0.010 0.001 

JES 0.003 0.007 

PDF 0.005 0.005 

color reconnection 0.001 0.004 

LO MC generator 0.005 0.005 



total 



0.024 0.017 



TABLE VI: Summary of inclusive asymmetries. Uncertainties 
include statistical, systematic, and theoretical uncertainties. 



sample 



level 



A' 



A^ 



data data 0.057 ± 0.028 0.073 ± 0.028 

mcQnlo tt+bkg 0.017 ± 0.004 0.001 ± 0.003 

data signal 0.075 ± 0.037 0.110 ± 0.039 

MC@NLO tt 0.024 ±0.005 0.018 ±0.005 

data parton 0.158 ± 0.074 0.150 ±0.055 

MCFM parton 0.058 ± 0.009 0.038 ± 0.006 



It is possible that the corrections in the presence of 
a large asymmetry would differ from the corrections de- 
rived from the symmetric PYTHIA. We have studied this 
possibility by applying the pythia- based response cor- 
rections to the OctetA model, which has an asymmetry 
like the data and a resemblance to the data in all other 
respects. We find that the bias in the corrected inclusive 
asymmetries is small, roughly 0.02, and we take this as 
evidence that the technique is essentially robust against 
perturbations of this kind. Since we have no reason to 
prefer the prediction of this or any other model, we do 
not include a modeling uncertainty. Our inclusive results 
assume the corrections and uncertainties calculated with 
the standard pythia model. 

Fig. [H] shows the Ay and yf^ distributions at all of the 
correction levels in the four-bin representation. The ef- 
fect of the background subtraction is clear. The tt signal 
(squares) derived from the background subtracted data 
can be directly compared with the pythia signal pre- 
diction, and continues to show the asymmetries. The 
corrected distribution at the parton-level (triangles) can 
also be compared to the symmetric pythia prediction. 



Table VI summarizes the measured asymmetries for 
the different levels of correction. It is interesting that 
at the data-level in the laboratory frame we compare to 
a model prediction that is consistent with zero. When 
the backgrounds are subtracted from the reconstructed 
data we can calculate the asymmetry for a pure ti sam- 
ple at the signal level, and compare directly to MC@NLO 
tt. The signal uncertainty here includes the uncertainty 
on the background correction. Correcting for acceptance 
and reconstruction resolution yields the tt parton-level 
asymmetry, where the uncertainty includes the effects 
listed in Table |V] The parton-level asymmetry may be 
directly compared with the standard model prediction of 

MCFM. 

The experimentally simple laboratory frame asymme- 
try exceeds the prediction by more than two standard 
deviations at all correction levels. The ti frame asymme- 
tries are similar in magnitude to the laboratory frame, 
but less significant because of the larger uncertainties. 
The ratio of the parton-level asymmetries in the two 
frames is A^^/A*^^ = 0.95 ± 0.41, where the error is cor- 
rected for the expected correlation across frames in the 
NLO QCD assumption. This measured ratio is consis- 
tent with the expected SM NLO value of 0.6, but the 
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uncertainty is large. 



B. Cross-Checks of the Inclusive Asymmetry 



Table [VTTI shows the asymmetries in the data when the 
sample is separated according to the Icpton flavor and the 
number of 6-tagged jets in the event. All of our simulated 
models predict asymmetries that are independent of the 
lepton type. Within the large errors, the data are con- 
sistent with this expectation. 

The 6-tagged sample contains 281 events with two b- 
tags. This double-tag sample is small, but has mini- 
mal backgrounds and robust jet-parton assignment. The 
double-tag sample is a special category of tt decays where 
both the b and b jet have | rj \< 1.0, but all of our simu- 
lation models predict similar asymmetries in single tags 
and double-tags. In the data the results are consistent 
across single and double-tags, albeit with reduced agree- 
ment in ^PP. We will discuss the double-tag consistency 
in the laboratory frame in more detail in Sec. |VIII E| 



TABLE VII: Measured asymmetries at the data-level for dif- 
ferent lepton and b-tag selections. 



selection 


A" 




inclusive 


0.057 


± 


0.028 


0.073 ±0.028 


electrons 


0.026 


± 


0.037 


0.053 ±0.037 


muons 


0.105 


± 


0.043 


0.099 ±0.043 


single fo-tags 


0.058 


± 


0.031 


0.095 ± 0.032 


double &-tags 


0.053 


± 


0.059 


-0.004 ±0.060 



VI. RAPIDITY DEPENDENCE OF THE 
ASYMMETRY IN THE tt REST FRAME 

In Sec. |IV| we discussed the importance of measur- 
ing the rapidity and M^j dependence of the asymme- 
try. The correlated dependence on both variables would 
be most powerful, but, given the modest statistical pre- 
cision of our current dataset, we begin with separate 
measurements of each. In this section we show how a 
Ay-dependence may be calculated from the results of 
Sec. VA The Mjf-dependence (as well as the correla- 
tion of and Ay) will be discussed in the sections 
following. 

In the standard model at NLO the tt frame asymme- 
try increases linearly with Ay, as seen in Fig. |6] The 
slope is significant, with the asymmetry reaching values 
of roughly 20% at large Ay. 

The Ay dependence of the asymmetry in our binned 
data can be calculated in each bin i of positive Ay as 



^"(Ay, 



jV(Ay,) - iV(-Ay,) 
iV(Ay,) ± iV(-Ay,) 



(6) 
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FIG. 6: Ay-dependence of A" according to mcfm. 
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FIG. 7: Parton level asymmetries at small and large Ay com- 
pared to SM prediction of mcfm. The shaded bands represent 
the total uncertainty in each bin. The negative going uncer- 
tainty for Ay < 1.0 is suppressed. 



A parton-level measurement of A*'(Ayi) in two bins 
of high and low Ay is available from the corrected Ay 
distribution in Fig. [5] We calculate the asymmetry sep- 
arately for the low rapidity difference inner bin pair 
I Ay I < 1.0 and the large rapidity difference outer bin pair 
|Ay| > 1.0. The systematic uncertainties in the bin-by- 
bin comparison are evaluated using the same techniques 
as in the inclusive measurement. Uncertainty in the back- 
ground shape and normalization assumptions cause a sig- 
nificant systematic uncertainty in the high Ay bin. 

The Ay-dependent asymmetries are shown in Table 
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TABLE VIII: The tt frame asymmetry A" at small and large 
rapidity difference, compared to the SM prediction of mcfm. 

sample level \Ay\ < 1.0 \Ay\ > 1.0 

data data 0.021 ± 0.031 0.208 ± 0.062 

data parton 0.026 ± 0.104 ± 0.056 0.611 ± 0.210 ± 0.147 
MCFM parton 0.039 ± 0.006 0.123 ±0.018 



|VIII| For the parton-level data, the first uncertainty is 
statistical and the second is systematic. The uncertainty 
on the MCFM prediction is dominated by the NLO theory 
uncertainty. For \Ay\ < 1.0, the small data-level asym- 
metry maps into a small parton-level value with large 
error. In the large Ay region the parton-level asymme- 
try is A"(|Ay| > 1.0) = 0.611 ± 0.270 (statistical and 
systematic errors added in quadrature) compared to the 
MCFM prediction of 0.123 ± 0.018. Fig. [7] displays the 
parton level comparison of asymmetries in data in the 
two Ay regions. 

VII. MASS DEPENDENCE OF THE 
ASYMMETRY IN THE tt REST FRAME 




300 400 500 600 700 800 900 

M„ {GeVIc') 



FIG. 9: Event distribution as a function of the total invariant 
mass Mu. 




M^-(GeV/c^; 



800 



FIG. 8: Mtt-dependence of A" according to MCFM. 



We now turn to the dependence of the asymmetry 
on the tt invariant mass M^j. The NLO QCD asym- 
metry also has a strong M^f dependence, as shown in 
Fig. [8j We generally expect the M^t dependence to con- 
tain characteristic information on the fundamental asym- 
metry mechanism. 

In this analysis, the value of M^j is derived from the 
same reconstruction used to compute the top quark ra- 
pidities. The Mjf distribution in our sample, shown in 
Fig. [9| is agreement with the standard prediction. Other 
recent studies of the top pair mass spectrum, including 



the parton-level differential cross section da/dMf^, show 
good agreement with the standard model [101 ISOl El] . 

Since the mass dependent behavior is usually described 
in the tt rest frame we focus on the asymmetry in rapidity 
difference Ay as a function of Mjj . The laboratory frame 
asymmetry derived with y^ is discussed in Sec. |VIII[ 

The underlying 2-dimensional distribution of Ay 
vs.Mff is shown on the left in Fig. [TO] We expect 
these variables to obey the simple kinematic relation- 
ship Mjj = 2mT cosh{Ay) , where ttit is the transverse 
mass of the tt system, and we see this in both the data 
and the prediction. It is clear that the prior measure- 
ment at large Ay captures only part of the region at 
large Mjj. Consequently, the separate measurements of 
the Ay- and Mtf-dependence of the asymmetry provide 
complementary information. 

Because cosh(Ay) is symmetric, this kinematic correla- 
tion is independent of the Mj^-dependence of any asym- 
metry in Ay. Because of the independence of my, the 
measurement at large |Ay| > 1.0 captures only part of 
the region at large Mjf. The separate measurements 
therefore provide complementary information. 

A mass dependent asymmetry ^'*(M(j j) is found by 
dividing the Ay — Mu plane into bins of mass Mjf ^ and 
calculating the asymmetry in each: 

,tt(^ . ^ NjAy > 0, Mu,, - N{Ay < 0, Mu.,) , . 
^ "^^^ N{Ay > 0, Mtu) + N{Ay < 0, M„- ,) ^ ' 

We use 50 GeV/c^ bins of Mtt below 600 GeV/c^, and 
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FIG. 10: Left: The Ay — Mu plane. Each dot represents one event, while the intensity of the shading shows approximately 
the event probability in the standard pythia based prediction. Right: The ti frame asymmetry in the data in bins of invariant 
mass Mtt, compared to the prediction of mc@nlo tt + backgrounds. The last bin includes all events with Mu > 700GeV/c^. 
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FIG. 11: The tt frame asymmetries in bins of invariant mass 
Mtt when the data is partitioned by lepton charge. 



100 GcV/c^ bins above that. The Mjj-dependent asym- 
metry in Ay is shown on the right in Fig. [lO] and Ta- 
ble |IXj compared to the prediction of mc@nlo in combi- 
nation with the standard background. The uncertainties 
in the plot are the statistical errors only; in the table the 
MC@NLO uncertainty contains both the statistical and 
theoretical component. In the bulk of the data at low 
mass the asymmetry is consistent with zero, while at high 
mass the asymmetry is consistently above the prediction. 
Fig.fTTIshows that when the data are separated by lepton 
charge, the asymmetries in the two independent samples 
behave in approximately opposite fashion. 



TABLE IX: The data-level asymmetry A" in bins of Mtt com- 
pared to the prediction of mc@nlo -I- backgrounds. 



bin-center 




A" 


(GeV/c^) 


N events 


data 


MC@NLO 


375 


532 


-0.019 ± 0.043 0.003 ±0.006 


425 


322 


-0.012 ± 0.056 


0.026 ±0.008 


475 


190 


0.158 ± 0.072 


0.013 ±0.010 


525 


95 


0.305 ± 0.097 


0.019 ±0.013 


575 


58 


0.138 ± 0.130 


0.063 ±0.020 


650 


34 


0.471 ± 0.151 


0.051 ±0.020 


750 


29 


0.103 ± 0.185 


0.091 ±0.022 



A. Asymmetries at High and Low Mass 



The large statistical errors in the ^*'(Mjf i) distribu- 



tion of Fig. 



10 do not allow any conclusion on the func- 
In order to make a quantitative mea- 



tional dependence, 
surement of A** (Mtt) in a simple, statistically meaningful 
way, we use a compact representation of A**(Mtj_j) into 
just two bins, below and above a given mass bound- 
ary. 

The boundary between the low and high mass regions 
is chosen based on a study of the color-octet samples de- 
scribed in the Appendix. These samples have yl"(Mtf i) 
distributions that are comparable to the data and reason- 
able for modeling the sensitivity in that variable. We find 
that the significance of the asymmetry at high mass is 
maximized when the bin division is at Mu = 450 GeV/c^, 
and therefore adopt this cut. 

Fig. [12] shows the Ay distributions when the data 
is divided into two regions, below and above Mu — 
450 GeV/c^. At low mass the asymmetry is consis- 
tent with zero. At high mass, the rapidity difference is 
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FIG. 12: Top: The distribution of Ay at low mass (left) and high mass (right). 
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FIG. 13: The distribution of Ay at high mass for events with negative leptons (left) and positive leptons (right). 



TABLE X: Charged and total asymmetries at the data-level, for all, low, and high M^f. 



selection 


all Mu 


Mti < 450 GeV/c" 


Mti > 450 GeV/c^ 


reco data 


0.057±0.028 


-0.016±0.034 


0.210±0.049 


MC@NLO 


0.017±0.004 


0.012±0.006 


0.030±0.007 




0.067±0.040 


-0.013±0.050 


0.210±0.066 




-0.048±0.039 


0.020±0.047 


-0.210±0.071 



broader, as expected from the kinematics, and an asym- 
metry is apparent. The top two lines of Table pC| compare 
the high and low mass asymmetries with the mc@nlo 
prediction. The uncertainty on the prediction combines 
the statistical and the theoretical uncertainties. At high 
mass the reconstructed asymmetry A** = 0.210 ± 0.049 
(stat) is more than three standard deviations above the 



prediction. 

The high mass Ay distributions for the two separate 
lepton charges are shown in Fig. [13] and the asymmetries 
in those distributions are summarized in the bottom part 
of Table |Xj Under the interchange of lepton charge, or, 
equivalently, under the interchange of t and f, the asym- 
metry at high mass is approximately reversed. This is 
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TABLE XI: Tests of the combined mass and rapidity correction procedure. True, reconstructed, and fully corrected asymmetries 
as found in the two mass regions. Uncertainties on predictions are statistical errors in the MC samples; at truth level these are 
negligible. 



Sample 


A" level Mtt < 450 GeV/c^ Mu > 450 GeV/c^ 


PYTHIA 


MC truth 0.002 


0.001 




reconstructed -0.011 ± 0.006 


-0.013 ±0.008 




corrected 0.001 ± 0.018 


0.006 ±0.014 


mcQnlo MC truth 0.043 


0.070 




reconstructed 0.015 ± 0.006 


0.043 ± 0.009 




corrected 0.066 ± 0.014 


0.086 ±0.011 


Octet A 


MC truth 0.081 


0.276 




reconstructed 0.024 ± 0.035 


0.183 ±0.010 




corrected 0.054 ± 0.022 


0.308 ±0.016 


Octet B 


MC truth 0.150 


0.466 




reconstructed 0.078 ± 0.036 


0.310 ±0.009 




corrected 0.187 ±0.024 


0.476 ±0.015 



consistent with CP conservation, and also a strong argu- 
ment against a false asymmetry arising in event selection 
or tt reconstruction, as neither the event selection nor 
reconstruction make reference to the lepton charge. 

The results here suggest that the modest inclusive 
asymmetry in the tt rest frame originates with a large 
asymmetry in a small population at high M^f. 



The Mass Dependent Asymmetry at the 
Parton-Level 



In the measurement of the inclusive asymmetry we 
used a simple matrix technique to correct the rapidity 
distributions for acceptance and resolution and derive 
parton-level asymmetries that could be compared with 
theory. We do this now for the mass dependent asymme- 
try in the ti frame. We divide the data into two bins in 
Ay, forward and backward, and two bins in mass, above 
and below 450 GeV/c^ and re- apply the well tested 4x4 
unfold machinery of the inclusive analysis. The proce- 
dure yields fully corrected, model-independent asymme- 
tries that can be compared with theoretical predictions. 

We represent the four bins of the parton-level dis- 
tribution of Ay and M^j by a single vector n = 
[nLF,nLB,nHF,nHB] where, for example, ulp is the 
number of forward events at low mass. As in the in- 
clusive case, we know that the true n distribution is 
modified by matrices representing the acceptance and 
then by the smearing in the reconstruction, so that 
?^signai = SAnparton- To measurc the parton-level value, 
we subtract backgrounds to recover the signal from the 
data, and then invert the transformation as in Eq. (5). 

As before, the matrices A and S are derived from 
Pythia Monte Carlo samples by comparing truth distri- 
butions to the same distributions after reconstruction. 
The bin-to-bin migration measured in the smearing ma- 
trix now includes the cross-terms between high and low 
mass and forward and backward Ay. The most signifi- 



cant migration is caused by mis-reconstructions that un- 
derestimate Mjf and smear the shape of the Mjf spectrum 
towards lower masses. 



TABLE XII: Systematic asymmetry uncertainties in the two- 
mass bin unfold 



Source 


Mtt < 450 GeV/c" Mtt 


> 450 GeV/c^ 


background size 


0.017 


0.032 


background shape 


0.003 


0.003 


JES 


0.005 


0.012 


ISR/FSR 


0.012 


0.008 


color reconnection 


0.009 


0.004 


PDF 


0.018 


0.004 


physics model 


0.035 


0.035 


total 


0.047 


0.049 



The accuracy of the procedure is first tested against 
simulated control samples using pythia and mc@nlo. 
The PYTHIA test uses a tt sample that is independent of 
the one used to create the response matrices. The top 



part of Table XI shows that the correction procedure is 
unbiased when operating on the symmetric pythia in- 
put. The MC@NLO sample allows us to study the accu- 
racy of the correction in measuring the NLO QCD effect. 
A small possible bias of ^ 0.02 is insignificant compared 
to the statistical uncertainty in the present data set. 

Next, we use the color-octet samples to test how well 
the correction derived from symmetric pythia can re- 
cover large parton-level asymmetries. The bottom half 
of Table |XI| shows that the correction procedure recov- 
ers both the high and low mass asymmetries to within 
a few percent of the true values. The corrections in the 
Octet sample show a possible ~ 0.02 — 0.03 bias that is 
marginally significant compared to the statistical preci- 
sion of the test. Because the Octet samples match the 
data well in the two key distributions Ay and Mtt (see 
Appendix) we expect that this is a representative mea- 
sure of possible model dependence in the correction, and 



17 



we assign a systematic uncertainty of 0.035 for this effect. 

Additional systematic uncertainties are evaluated in 
a manner similar to the inclusive case. These uncertain- 
ties are estimated by repeating the analysis while varying 
the model assumptions within their known uncertainties 
for background normalization and shape, the amount of 
initial- and final-state radiation (ISR/FSR) in pythia, 
the calorimeter jet energy scale (JES), the model of fi- 
nal state color connection, and parton distribution func- 
tions (PDF). Table XII shows the expected size of all 
systematic uncertainties. The physics model dependence 
dominates. 



TABLE XIII: 
to prediction. 



Asymmetry yl" at high and low mass compared 



selection 


Mtt < 450 GeV/c" 


Mti > 450 GeV/c^ 


data 


-0.016 ±0.034 


0.210 ±0.049 


ti+hkg 


+0.012 ±0.006 


0.030 ±0.007 


(mcQnlo) 






data signal 


-0.022 ±0.039 ±0.017 


0.266 ± 0.053 ±0.032 


ti 


±0.015 ±0.006 


0.043 ± 0.009 


(mc@nlo) 






data parton 


-0.116 ±0.146 ±0.047 


0.475 ±0.101 ±0.049 


MCFM 


±0.040 ± 0.006 


0.088 ±0.013 
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FIG. 14: Parton-level asymmetry in Ay at liigh and low mass 
compared to mcfm prediction. The shaded region represents 
the total uncertainty in each bin. 



Table XIII compares the low and high mass asymme- 
try to predictions for the data level, the background sub- 
tracted signal-level, and the fully corrected parton-level. 
The MC predictions include the 15% theoretical uncer- 
tainty. At low mass, within uncertainties, the asymmetry 
at all correction levels agrees with predictions consistent 
with zero. At high mass, combining statistical and sys- 
tematic uncertainties in quadrature, the asymmetries at 
all levels exceed the predictions by more than three stan- 
dard deviations. The parton-level comparison is summa- 
rized in Fig.nS For Mtt > 450 GeV/c^, the parton-level 
asymmetry at in the ti rest frame is A" = 0.475 ±0.114 
(stat±sys), compared with the MCFM prediction of 
A" = 0.088 ±0.013. 
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FIG. 15: Distribution of reconstruction x^- Black crosses 
are data, histogram is sig±bkg prediction. The last bin on the 
right contains all events with > 100. 



VIII. CROSS-CHECKS OF THE MASS 
DEPENDENT ASYMMETRY 

The large and unexpected asymmetry at high mass de- 
mands a broader study of related effects in the tt data. 
We look for anomalies that could be evidence of a false 
positive, along with correlations that could reveal more 
about a true positive. In order to avoid any assumptions 
related to the background subtraction, we make compar- 
isons at the data level, appealing when necessary to the 
full ti + bkg simulation models. 



A. Lepton Type 

All of our simulated models predict asymmetries that 
are independent of the lepton type: pythia predicts 
asymmetries that are consistent with zero, and the Octet 
models predict asymmetries that are consistent with each 
other. The data are shown in Table |XlV] At high mass, 
both lepton types show positive asymmetries consistent 
within errors. 
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TABLE XIV: Data level asymmetries A" for different event selections. In the case of no-b-fit, the tt reconstruction has been 
run without the constraint that b-tagged jets be associated with b-partons. 



selection 


N events 


all Mtt 


Mtt < 450 GeV/c" M 


i > 450 GeV/c" 


standard 


1260 


0.057±0.028 


-0.016±0.034 


0.212±0.049 


electrons 


735 


0.026±0.037 


-0.020±0.045 


0.120±0.063 


muons 


525 


0.1U5±0.U4.3 


-U.0i2±0.0o4 


0.348±0.080 


data < 3.0 


338 


0.030±0.054 


-0.033 ±0.065 


0.180 ±0.099 


data no-b-fit 


1260 


0.062±0.028 


0.006 ± 0.034 


0.190 ±0.050 


data single b-tag 


979 


0.058±0.031 


-0.015±0.038 


0.224±0.056 


data double b-tag 


281 


0.053±0.059 


-0.023±0.076 


0.178±0.095 


data anti-tag 


3019 


0.033±0.018 


0.029±0.021 


0.044±0.035 


pred anti-tag 




0.010±0.007 


0.013±0.008 


0.001±0.014 


pre-tag 


4279 


0.040±0.015 


0.017±0.018 


0.100±0.029 


pre-tag no-b-fit 


4279 


0.042±0.015 


0.023±0.018 


0.092±0.029 



B. Reconstruction 

It is conceivable that a reconstruction error could pro- 
duce an asymmetry from symmetric inputs. The qual- 
ity of the reconstruction is summarized by a that 
measures the consistency of the solution with the tt hy- 
pothesis. The distribution of in our sample, shown in 



Fig. 15 is in very good agreement with the prediction, 
including a good match on the long tail. When the sam- 
ple is restricted to high quality fits with x^ < 3.0, we find 
338 events in which A" = —0.033 ± 0.065 at low mass 
and A" = 0.180 ± 0.099 at high mass. Although the 
statistical precision is diminished in this small sample, 
it suggests that the high mass asymmetry is present in 
the best reconstructed events. Since the x^ requirement 
rejects a significant fraction of the background, it also 
suggests that the high mass asymmetry is not a back- 
ground related effect. 

To test for possible reconstruction biases related to b- 
tagging, we re-run the reconstruction algorithm removing 
the constraint that 6-tag jets be matched to b partons. 
We find ^" = 0.006 ± 0.034 at low mass and A" = 
0.190 ± 0.050 at high mass. When we further separate 
the events by lepton charge, the Ayih asymmetries are 
A" = -0.190 ± 0.074 and A+ = 0.190 ± 0.069. The 



Ih 



large forward-backward charge asymmetry at high mass 
is seen to be independent of the use of b-jet identification 
in the reconstruction. 



In the background dominated anti-tags, the inclusive 
and low mass samples have small asymmetries that agree 
with the prediction. In the high mass anti-tag sample we 
find A*-*- = 0.044± 0.035, consistent with either the model 
prediction of zero or a slight excess due to the tt compo- 
nent there. Mixing backgrounds and tt in the expected 
ratio and assuming the tt component has an asymme- 
try of 0.266 (as in Table XIII), we find a total expected 
asymmetry in the anti-tag sample of A" — 0.079 ± 0.034 
in agreement with the data. 

The lepton+jets sample with no 6-tagging is the "pre- 
tag" sample. Our standard pythia + background model 
predicts pre-tag asymmetries consistent with zero for all 
mass categories. The asymmetries in the data are shown 



in Table XIV At low mass the asymmetry in the pre-tags 
is consistent with zero. At high mass, the pre-tag sample 
has a significant asymmetry 0.100 ± 0.029. If we assume 
that tt signal at high mass has = 0.266 as in Ta- 



ble XIII and combine tt with our standard backgrounds 



in the expected pre-tag ratio, we predict a pre-tag asym- 
metry of v4** ~ 0.111 ± 0.028, in good agreement with 
the data. 

As a final check in the pre-tag sample, we repeat the 
exercise of running the reconstruction without the con- 
straint that 6-tagged jets are used as b-partons. The 
results are shown in the bottom row of Table IXIVI The 
asymmetry at high mass is 0.092 ± 0.029, a significant 
effect in a sample that makes absolutely no reference to 
^-tagging. 



C. B-Jet Identification 

All of our simulated models predict asymmetries that 
are independent of whether one or two jets are 6-tagged. 
In single and double 6-tagged samples pythia predicts 
asymmetries that are consistent with zero, and the Octet 
models predict asymmetries that are consistent with each 
other. In the data, the two cases are consistent with each 
other, although the statistical precision on the double 
tagged sample is marginal. 



D. Jet Multiplicity 



In Sec. |IV A| we discussed the two components of 
the NLO QCD asymmetry: (1) radiative corrections to 
quark-antiquark production and (2) interference between 
different amplitudes contributing to the ttj final state. 
The two contributions have opposite signs. At NLO, 
the first is positive and dominant for the inclusive mea- 
surement, while the second is negative and subdominant. 
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TABLE XV: mc@nlo predictions for A" in reconstructed tt signal (no 
backgrounds) as a function of A/jf and jet multiplicity. The uncertainties 
reflect MC statistics only. 

selection all Mtt < 450 GeV/c^ Mtt > 450 GeV/c^ 

inclusive 0.024 ± 0.004 0.015 ± 0.005 0.043 ± 0.007 

4- jet 0.048 ±0.005 0.033 ± 0.006 0.078 ± 0.009 

5- jet -0.035 ±0.007 -0.032 ± 0.009 -0.040 ±0.012 



TABLE XVI: Asymmetries A" in the data as a function of jet multiplicity. 

selection N events all Mtt < 450 GeV/c^ Mu > 450 GeV/c^ 

inclusive 1260 0.057±0.028 -0.016±0.034 0.212±0.049 

4- jet 939 0.065±0.033 -0.023±0.039 0.26±0.057 

5- jet 321 0.034±0.056 0.0049±0.07 0.086±0.093 



TABLE XVII: Reconstruction level asymmetries A^^ in the laboratory frame. 



selection 


all Mu 


Mtt < 450 GeV/c^ M 


„- > 450 GeV/c" 


data reco 


0.073±0.028 


0.059±0.034 


0.103±0.049 


MC@NLO 


0.001±0.003 


-0.008±0.005 


0.022±0.007 


At 


-0.070±0.040 


-0.028±0.050 


-0.148±0.066 




0.076±0.039 


0.085±0.047 


0.053±0.072 


single b-tags 


0.095±0.032 


0.079±0.034 


0.130±0.057 


double b-tags 


-0.004±0.060 


-0.023±0.076 


0.028±0.097 



Since only the second term produces ttj events, we ex- 
pect that the QCD asymmetry will be a function of the 
jet multiplicity. 

We have studied the jet multiplicity dependence of A** 
in MC@NLO. We define 4-jet events as those with four 
jets with Et > 20 GeV and | ?7 |< 2.0 and no other such 
jets. We define 5-jet events as those with at least five jets 
with Et > 20 GeV and | rj |< 2.0. The mcOnlo predic- 
tion for the pure tt signal after reconstruction is shown 
in Table |XV[ The 5-jet asymmetries are negative, as ex- 
pected. The exclusive 4-jet sample shows asymmetries 
that are roughly double those in the inclusive sample. 

As we discussed in Sec. |IV A} the reliability of the NLO 
picture has recently been called into question by NNLO 
calculations of the ttj component [25] , which reduce the 
negative asymmetry there to close to zero. However, 
since no NNLO calculation exists for the exclusive 4-jet, 
inclusive, or mass dependent asymmetries, the MC@NLO 
prediction in Table [XV| remains our comparison point. 

The jet multiplicity dependence of the asymmetries in 
the data is shown in Table |XVl| Vetoing events with ex- 
tra jets does not produce a significant increase in the 
asymmetry. However, in the 5-jet sample, the high mass 
asymmetry is consistent with zero. With a larger sample 
and better precision it might be possible to use the jet 
multiplicity to test whether the observed asymmetry is 
an amplified version of the QCD charge asymmetry or a 



different effect altogether. 



E. Frame Dependence 

As in the inclusive analysis, it is interesting to compare 
A" to APP. In the NLO QCD e ffect, t he frame depen- 
dence of the asymmetry (see Sec. IV A) persists at high 
mass. For M^f > 450 GeV/c^ our mc@nlo model pre- 
dicts the ratio of reconstructed asymmetries in the two 
frames App/A^^ ~ 0.74. The OctetA model predicts less 
mass dependence, with a ratio of 0.90. 

The lab frame data asymmetries above and below 
Mtt = 450 GeV/c^ are shown in Table 



XVII 



The varia- 
tion of the asymmetry across the 450 GeV / mass edge 
is not as distinct as in the tt frame, and the deviation 
from the mc@nlo prediction is not as significant. Within 
the large errors, the asymmetries in the two lepton charge 
samples are consiste nt with CP invariance. 

Comparing Tables [xVH] and [xj the ratio of App to A" 
at high mass is 0.49±0.21, lower than both the mc@nlo 
and Octet models. We have used pseudo-experiment 
techniques to evaluate the statistical consistency of this 
ratio with the models, using a large number of simulated 
experiments that differ by Poisson fluctuations in the Ay 
and —qyii distributions. A App/A" ratio of 0.49 or less 
occurs in 14% of pseudo-experiments with mc@nlo, but 
in < 1% of experiments with OctetA. 
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Finally, we look at A^^ as a fun ction of the 6-tag mul- 
tiplicity. We observed in Sec. VII that the inclusive A^^ 
is zero in the double 6-tagged events. In Table pCVlI[ we 
see that this pattern persists at high mass, although the 
statistical precision is poor. Appealing again to pseudo- 
experiments with Poisson fluctuations, we find that a ra- 
tio of double to single tag A^^ as small as that in the data 
occurs in 6% of all pseudo-experiments with MC@NLO. 
We conclude that the low value of A^p in the double b- 
tagged sample is consistent with a statistical fluctuation. 



IX. CONCLUSIONS 

We have studied the forward-backward asymmetry of 
top quark pairs produced in 1.96 TeV pp collisions at 
the Fermilab Tevatron. In a sample of 1260 events in 
the lepton-f-jet decay topology, we measure the parton- 
level inclusive asymmetry in both the laboratory and tt 
rest frame, and rapidity-dependent, and Mjf-dependent 
asymmetries in the tt rest frame. We compare to NLO 
predictions for the small charge asymmetry of QCD. 

The laboratory frame measurement uses the rapidity 
of the hadronically decaying top system and combines 
the two lepton charge samples under the assumption of 
CP conservation. This distribution shows a parton-level 
forward backward asymmetry in the laboratory frame of 
APP = 0.150 ± 0.055 (stat-hsys). This has less than 1% 
probability of representing a fluctuation from zero, and 
is two standard deviations above the predicted asymme- 
try from NLO QCD. We also study the frame-invariant 
difference of the rapidities. Ay — yt — yt, which is pro- 
portional to the top quark rapidity in the tt rest frame. 
Asymmetries in Ay are identical to those in the t pro- 
duction angle in the tt rest frame. We find a parton-level 
asymmetry of A" = 0.158 ± 0.075 (stat-t-sys) , which is 
somewhat higher than, but not inconsistent with, the 
NLO QCD expectation of 0.058 ± 0.009. 

In the tt rest frame we measure fully corrected asym- 
metries at small and large Ay 

Ay| < 1.0) ^ 0.026 ± 0.118 
A"(|Ay| > 1.0) = 0.611 ± 0.256 

to be compared with mcfm predictions of 0.039 ± 0.006 
and 0.123 ± 0.008 for these Ay regions respectively. 

In the tt rest frame the asymmetry is a rising function 
of the tt invariant mass Mjf, with parton level asymme- 
tries 

>l"(Mtf < 450 GeV/c2) = -0.116 ± 0.153 
{Mtt > 450 GeV/c^) = 0.475 ± 0.114 

to be compared with mcfm predictions of 0.040 ± 0.006 
and 0.088 ± 0.013 for these Mft regions respectively. The 
asymmetry at high mass is 3.4 standard deviations above 
the NLO prediction for the charge asymmetry of QCD, 
however we are aware that the accuracy of these theo- 
retical predictions are under study. The separate results 



at high mass and large Ay contain partially independent 
information on the asymmetry mechanism. 

The asymmetries reverse sign under interchange of lep- 
ton charge in a manner consistent with CP conservation. 
The tt frame asymmetry for Mjf > 450 GeV/c^ is found 
to be robust against variations in tt reconstruction qual- 
ity and secondary vertex b-tagging. When the high-mass 
data is divided by the lepton flavor, the asymmetries 
are larger in muonic events, but statistically compatible 
across species. Simple studies of the jet multiplicity and 
frame dependence of the asymmetry at high mass may 
offer the possibility of discriminating between the NLO 
QCD effect and other models for the asymmetry, but the 
statistical power of these comparisons is currently insuf- 
ficient for any conclusion. 

The measurements presented here suggest that the 
modest inclusive tt production asymmetry originates 
from a significant effect at large rapidity difference Ay 
and total invariant mass Mtt- The predominantly qq 
collisions of the Fermilab Tevatron are an ideal environ- 
ment for further examination of this effect, and additional 
studies are in progress. 
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X. APPENDIX: THE COLOR-OCTET MODELS 

In the generic color-octet model of Rcf. |8] , the gluon- 
octet interference produces an asymmetric cos(0*) term 
in the production cross section. The couplings of the 
top and the light quarks to the massive gluon have op- 
posite sign, giving a positive asymmetry as seen in the 
data. This was implemented in the madgraph frame- 
work, and the couplings and Mq were tuned to reason- 
ably reproduce the asymmetries and Mjj distribution of 
the data [26] . The sample called Octet A, with couplings 
9v = 0, gA{q) = 3/2, gA{t) = -3/2, and mass Ma = 2.0 
TeV/c^, has parton level asymmetries of A^p = 0.110 and 
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FIG. 16: The distributions of Mu (left) and Ay (right) in the OctetA sample compared to pythia. 
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FIG. 17: The reconstructed M dependence of yl** in the color-octet models. Left: OctetA. Right: OctetB. 



TABLE XVIII: A" and significance (in units of standard deviation) in Octet samples for events with Mtt above the bin-edge. 
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A*"*" — 0.157. The LO cross section for this sample is 6.1 
pb, in good agreement with the LO madgraph cross 
section for standard model production at 6.0 pb. 

The sample is showered with pythia, run through the 
CDF-II detector simulation, and then subjected to our 
selection and reconstruction. The M^t and Ay distribu- 
tions in OctetA are compared to the pythia versions 
in Fig. [16] The M^j distribution is a good match to 
PYTHIA, and we have checked that related transverse 
variables are also well-modeled. The rapidity distribu- 
tions after selection and reconstruction have asymmetries 
APP = 0.073 ± 0.006, A" = 0.079 ± 0.006, which are rea- 



sonable matches to the data. 

The complementary OctetB sample has the same cou- 
plings but Mc = 1.8 TeV/c^, giving parton-level asym- 
metries APP = 0.205 and ^" = 0.282. The tt cross sec- 
tion increases by 5% and the reconstructed mass distri- 
bution has a slight excess at the high mass relative to 
PYTHIA. The data level asymmetry A" = 0.16 ± 0.006, 
significantly higher than the data. 

These models both show the same approximate factor 
of 2 dilution between data-level and parton-level asym- 
metries that is seen in the data and in mc@nlo. 

Since these models have relatively low-lying octet 
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FIG. 18: Top: The distribution of Ay at low mass (left) and high mass (right). The pythia baseline model is the filled 
histogram and the Octet samples are the points. 



TABLE XIX: Reconstructed asymmetry below and above = 450 GeV/c^ in Octet models. 

sample aU Mtt < 450 GeV/c^ Mtt > 450 GeV/c^ 

OctetA 0.080 ±0.006 0.024 ±0.007 0.180 ±0.010 

Octets 0.160 ±0.006 0.078 ± 0.007 0.310 ± 0.009 



masses near 2 TeV/c^ we expect a significant Mii- 
dependent asymmetry over our experimental range. The 
A**(Mjt_j) behavior for the two color-octet samples is 
shown in Fig. |17[ Both show a smooth and significant 
rise of the asymmetry with increasing mass. 

In Sec. |VII A| we discussed a simple representation of 
A** (Mjf ) with two regions of low and high M^f. The ques- 
tion for that representation is how to choose the bound- 
ary mass between high and low. Table XVIII shows 
the asymmetry, uncertainty, and significance ^"/tr^tt at 
high mass as a function of the mass threshold for both 
octet models. The uncertainties are calculated assuming 
the data sample size of 5.3 fb^"'^. For both samples, the 
significance of the asymmetry at high mass is maximum 
at reconstructed Af^f = 450 GeV/c^. Looking at Fig. 9j 
we see that this is reasonable: 450 GeV/c^ cuts oflf t' 



bulk of the low mass peak while retaining good statistics 
on the tail. 

Fig. [18] compares the Ay distributions in the OctetA 
sample and pythia when the events are divided into sam- 
ples below and above Mfi = 450 GeV /<? . The Ay distri- 
bution is much broader at high mass, as expected, and 
shows a marked asymmetry. 

The reconstructed asymmetries at high and low mass 
in the color octet samples are given in Table |XIX[ The 
uncertainties here reflect the Monte Carlo statistics. At 
high mass the color octet samples have large asymmetries 



as seen in Fig. 18 At low mass, the models have small but 
significant asymmetries, especially in OctetB. We have 
checked that these asymmetries are charge asymmetries, 
reversing sign under interchange of lepton charge. 
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